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Abstract 



Three binaries are now established sources of emission at very high 
energies (>10 11 eV). They are composed of a massive star and a com- 
■ pact object. The emission can be due to the interaction of the rel- 

ativistic wi nd from a young ms puls ar with the stellar wind of the 
companion ()Maraschi fc Trevest H98l]b by which rotation-power ends 
up as non-thermal flux. Variations at VHE energies are explained as 
due to 77 absorption and/or changes in shock location along the orbit. 
Resolved radio emission is due to cooling particles trailing the pulsar. 

o' 

1 The observational backdrop 

Thanks to the new generation of Cherenkov telescopes HESS and MAGIC, 
the past two years have seen a spectacular increase in the number of 
known sources emitting photons with energies greater than a TeV (see 
the contributions by Cui, Santangelo and Teshima in this volume). There 
are now about 40 TeV sources, mostly concentrated in the Galactic 
plane. Three of the Galactic sources have been associated with binary 
stars on the b asis of their positional c oincidence, variability and spectra : 

" J2005hh . LS 5039 (jAharonian et al.l .l2005al 



PSR B1259-63 ( Aharonian ct al 



and LS I +61°303 (jAlbert et al 



, 2006). 
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1.1 The binaries 



Those three gamma-ray binaries are composed of a massive star (O or Be 
type) and a compact object, neutron star or black hole, in an eccentric orbit. 
The orbital periods are ps 4 days (LS 5039), 26 days (LS I +61°303) and 
1237 days (PSR B1259-63). Given the small mass functions, radial velocity 
measurements have been unable to pin down the nature of the compact 
object in LS 5039 and LS I +61°303. In both cases, a black hole would imply 
seeing the systems at i ~ 30°, while a neutron star would imply a more edge- 
on i ?s 60°. Regardless, the presence of a luminous early- type star, providing 
a copious source of seed photons for inverse Compton scattering, is probably 
instrumental to very high energy (VHE) 7-ray emission. 

With distances ranging from 1.5-3 kpc, the TeV emission is at a level 
~ iq33-34 er g g — 1^ comparable to the X-ray emission from these sources. 
Remarkably, all have been detected in the radio, a feature otherwise shared 
by only a handful of high-mass X-ray binaries. Bar from periodic radio out- 
bursts in LS I +61°303 and PSR B1259-63, very likely related to passage of 
the compact object through the dense equatorial wind of the Be companion 
(LS 5039 has a 09V companion - and no such outburst), the emission from 
the systems is surprisingly steady up to timescales of years. This contrasts 
with typical HMXBs. Their spectral energy distributions superpose well 
with each other, with a rising spectral luminosity from radio to 0.1-1 MeV, 
a flat spectrum up to 1 — 10 GeV and a drop at TeV energies. 

1.2 What powers the VHE emission? 

The high-energy emission could be powered by accretion via capture of the 
intense stellar wind from the companion, but the Bondi-Hoyle rates are low 
and observational signatures of accretion are lacking. Still, microquasar-type 
models have been proposed in which the 7-rays are emitted in a relativistic 
jet (see contribution by Paredes in this volume). Resolved compact radio 
emission in LS 5039 and LS I +61°303 on scales of 10-100 mas and inter- 
preted as jet emission, constitutes a strong motivation. 

On the other hand, the spectral and temporal similarities between the 
three binaries detected in VHE 7-rays hint at a common scenario. This is 
naturally provided for by PSR B1259-63. In this system, the relativistic wind 
emitted by the young 48 ms pulsar is contained by the stellar wind. Particle 
acceleration at the termination shock leads to X-ray synchrotron emission 
and 7-ray inverse Compton emission with the star photons. The measured 
pulsar spindown power is a modest 8 • 10 35 erg s _1 , implying rotation-power 



is converted into radiation with a high efficiency. 

The X-ray spectra have photon index T ~ 1.5, harder than those of Crab- 
like high-luminosity pulsars but comparable to those of pu lsar wind nebul a 
(PWN) powered by young pulsars with E ~ 10 36 erg s -1 (|Gottheld . 12003). 
X-ray pulses might be detectable underneath the PWN emission if scattering 
is negligible. However, pulsed radio emission would be strongly absorbed by 
the dense stellar wind in LS 5039 and LS I +61°303, with free-free opacities 
r > 10 3 . The pulsed radio emission does vanish in PSR B1259-63 around 
periastron, when the pulsar probes the densest region of the Be wind. At 
that point, the orbital separation is comparable to the maximum separation 
in LS I +61°303, and much larger than in the very compact LS 5039; 

The advantages of the PWN scenario, discussed below and in Dubusl 
(2006a), are that it readily explains the spectrum, level and stability of 
the emission without requiring additional assumptions from what is already 
known from plerions, and that it offers a very fruitful, common framework 
for interpretation. 



2 High energy emission from a compact PWN 

After the termination shock, a nebula of accelerated particles forms behind 
the pulsar, in the direction opposite to the apparent motion. Numerical 
simulations show the flow can be quite complex, but the major factor in 
setting the properties of the nebula is the distance of the termination shock 
Rs from the pulsar. R$ is set by the balance of the ram pressures of the 
stellar and pulsar winds. 



2.1 The model 

A straightforward approach to test a PWN scenario for LS 5039 and 
LS I +61°303 is to assume a pulsar with properties similar to those of 
PSR B1259-63 but to change the stellar wind and orbit according to the 
optical observations of the binaries. The pulsar wind is therefore assumed 
to carry an energy ~ 10 36 erg s , composed of mono-energetic e + e~ pairs 
with Lorentz factor 7^ ~ 10 5 and a small magnetic field so that the ratio 
of magnetic to kinetic energy density a ~ 10~ 2 , all reasonable values for a 
young, spinning-down pulsar. 

For LS 5039, the wind mass-loss rate is ~ 10 7 M yr _1 with a terminal 
velocity ~ 2000 km s _1 so that 

E/4ttR 2 s c = Pw v 2 w -> R s « 1.5 • 10 n do.i^36 2 cm i 1 ) 



where do.i = 0.1 AU is the orbital separation. The termination shock occurs 
close to the pulsar at a small fraction of the orbital separation. 

The principles of an emission model are easily derived by analogy with 
plerionic supernova remnants. MHD shock conditions set the magnetic field 
intensity, flow speed (« c/3), density etc. Particles are assumed to be ac- 
celerated to a dN oc 7 _2 c?7 power-law up to the energy at which radiative 
losses occur on a timescale smaller than Bohm diffusion. Density conser- 
vation then sets the minimum 7 of the distribution. The synchrotron and 
inverse Compton emission on the stellar photon field (which occurs in the 
Klein-Nishina regime) can then be calculated. Particles cool as they move 
away from the pulsar; changes in flow conditions are given by a Bernouilli 
equation. 



2.2 Spectral energy distribution 

The major features of the spectral energy distribution derive from sim- 
ple considerations on timescales. The acceleration timescale for TeV elec- 
trons (76=10 6 ) is t acc rs 0.06 7e/-Bi s with B\=\ G the typical magnetic 
field intensity. The synchrotron cooling timescale is t$ ~ 770/-B 2 76 s. 
Inverse Compton losses on stellar photons occur on a timescale tic ~ 
20 76<^ 1/ [In 7 6 + 1.4] (T* i4 .R*,io) s with T*=40,000 K and i^=10 R Q ap- 
propriate for the O star in LS 5039. 

At the highest energies, 7 max is set by synchrotron losses (iacc=^s)> which 
dominate over IC losses above a critical 7b r k given by (ts=iic) : 

7brk « 6 • 10 6 (T* i4 /2*,io)/(5ido.i). (2) 

Assuming continuous injection of electrons with 7 -2 spectrum, the steady- 
state distribution is steepened by radiative losses above 7t, r k to 7 -3 . This 
produces a flat synchrotron spectrum (in vF v ) above the frequency z^s 
corresponding to 7t, r k (Eq. 3). The inverse Compton spectrum in Klein- 
Nishina regime roughly reflects the particle distribution with vF v oc v 
( Moderski et al. . 20051 ) above /iz^ic=7brk m eC 2 (Eq. 4). Inefficient Klein- 



Nishina losses dominate below Tbrk, producing a hard injection-like spec- 
trum. Hence, the synchrotron spectrum below 1/3 is vF v oc i^ ' 5 and the 
inverse Compton spectrum below v\q is flat. The characteristic frequencies 
are given by 

hva « 750(r M i?*,ioM).i) 2 ABi keV (3) 
hv lc « 4(T M i2*,i /do.i)/Si TeV. (4) 




Figure 1: Spectral energy distributions for the 7-ray binaries with PWN emission 
models. For LS 5039, the pulsar wind has £=10 36 erg/s, j w =5 10 4 and cr=0.01. 
The termination shock is at i? s =4 10 11 cm at apastron (solid line) and 2 10 11 cm at 
periastron (dashed line). The plot includes VHE 77 absorption (Fig. 2) but does 
not include emission from an associated cascade. For LS I +61°303, E—10 36 erg/s, 
j w =5 10 4 , R s =2 10 11 cm and tr=0.02 at periastron in the dense equatorial wind 
(dashed line), changing to i? s =4 10 12 cm and cr=0.005 at apastron in the polar 
wind (solid line). The small R s at periastron implies a high B so enhanced X- 
ray synchrotron and lower IC at VHE energies. For PSR B1259-63, the periastron 
parameters (dashed line) are i? s =10 12 cm, E=8 10 35 erg/s, 7^=4 10 5 and er=0.005. 



These considerations are quite general and will apply to any leptonic 
model close to the compact object. Numerical SEDs for the three binaries 
are shown in Fig. 1, consistent with the above expectations. The radio to 
MeV fluxes are in good agreement with observations. EGRET fluxes are 
underestimated but uncertainties in the GeV Galactic background emission 
could explain the discrepancy (cascades too). The TeV is discussed below. 



2.3 Periodicity in LS 5039 

An attractive feature of the PWN model is that the evolution of the SED 
with orbital phase is very easy to understand. If the stellar wind does 
not vary along the orbit, the magnetic field intensity is tied to Rs so that 
B cx 1/ R s oc 1/d (Eq. 1). Hence, the spectrum should be unchanged at TeV 



energies but the hard X-ray spectral break should move as 1/d. 

In LS 5039, d varies only by a factor 2, so few changes are expected 
in the intrinsic SED. However, TeV photons embedded deep in the stellar 
radiation field can pair produce with the UV photons, leading to periodic 
absorption with maximum flux at inferior conjunction (4>=0. 7, Fig. 2), whe n 
the observed TeV photons travel directly away from the star ( DubusL 2006b). 
A few weeks after giving this talk, HESS reported that the TeV flux is indeed 
strongly modulated on P OT b, confirming this basic picture and exc luding 
emission far down (>1 AU) a relativistic jet ( Aharonian et al. . 2006) . The 



modulation is stable, as expected in the PWN model. 

Spectral changes are observed with the modulation (Fig. 1). At low 
fluxes (~ periastron), the spectrum is steep. The source is never fully 
absorbed: VHE absorption must initiate an electro-magnetic cascade re- 
emitting 7-rays at lower energies, notably in the EGRET band. Obvi- 
ously, this complicates the interpretation. At high fluxes (~ apastron where 
there should be little absorption), the TeV spectrum is flat with a break 
at « 8 TeV, suggesting 7b r k changes strongly with phase. The geometry of 
the interaction region may play a crucial role here. For instance, a value of 
S=0.25 G (instead of 4 G in Fig. 1) would produce a 8 TeV break but im- 
plies a very low-a pulsar wind reaching to the star surface. Alternatively, an 
additional component might produce this hard spectrum, perhaps hadronic: 
theoretical investigations show nuclei in the pulsar wind are necessary for 
significant non-thermal acceleration. 



2.4 Periodicity in LS I +61°303 and PSR B1259-63 

MAGIC recently reported detecting LS I +61°303 at apastron but did not 
detect it at phases </>=0-0.1 when the pulsar is in front of the star i.e. when 
absorption is actually weakest 1 (Fig. 2). In LS I +61°303, the wider orbit and 
orientation are such that the opacity from the pulsar towards the observer 
is at most 2, while it reaches 40 in LS 5039. 

The TeV variability in LS I +61°303 is more likely due to the presence 
of a Be star than to 77 absorption. Be stars have a slow, dense equatorial 
wind and a fast, polar wind (similar to that of the O star in LS 5039). The 
equatorial wind extends up to ~ 20 R* so that the pulsar in LS I +61°303 and 
PSR B1259-63, who both have Be companions, plunges periodically through 
the dense material. This crushes the PWN, like the Earth magnetosphere 
during a Solar storm. The termination shock is much closer to the pulsar 

1 Gupta & Bottcher (2006) find 77 absorption explains the non-detection (at 0=0-0.1) 
using 0=0.93, i.e. superior conjunction where pair production is maximal. 
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Figure 2: Expected modulation of the integrated VHE flux above 100 GeV due to 
pair production on stellar photons in LS 5039 (left) and LS I +61°303 (right). For 
LS 5039, peak emission is at 0=0.7, as observed by HESS. Non-zero TeV fluxes are 
also detected at 0=0-0. 2 (Fig. 1), probably due to cascading. For LS I +61°303, 
MAGIC detects the source at 0=0.3-0.5 but sees no flux at phases 0-0.1 despite 
little absorption. The variability is more likely to be intrinsic. (Note that periastron 
is often defined set at 0=0.23 instead of 0=0 in LS I +61°303.) 



so B is much higher, synchrotron losses are greater and the IC emission 
lower. If B is high, v\q is low and the steep IC spectrum leads to fainter 
TeV emission. At apastron, the stellar wind is tenuous, the termination 
shock further away from the pulsar, B is low, 7t> r k higher and IC emission 
more importan t at TeV, just as observed by MAGIC. This was predicted in 
Duhns! (j2006ah . 



In PSR B1259-63 the HESS lightcurve indicates a more complex, asy- 
metric behaviour around periastron. An important thing to note when com- 
paring the two systems is that the periastron separation in PSR B1259-63 
is close to the apastron separation in LS I +61°303, and that the Be disk 
appears inclined with respect to the orbital plane in PSR B1259-63. There- 
fore, PSR B1259-63 samples the equatorial disk twice around 0=0 and at 
comparatively lower densities; IC emission would first decrease and then re- 
cover after each passage. This works well: the X-ray and TeV emission show 
peaks after the disc crossings inferred from radio dispersion measures 2 . 

Further complications in computing the high energy spectra and 
lightcurves for the binaries may arise from IC anisotropy: the respective 



2 The peaks are less li kely to be due to the actual Be disc crossing, as proposed by 
IChernvakova et al.l (|200(t) . The required orientation of the Be disc with respect to the 
orbit would imply the 2 nd crossing occurs 30 days after periastron, inconsistent with the 
pulsar radio eclipse ceasing 10 days earlier than that. Furthermore, the compact object is 
then 40 i?* away from the star, beyond the truncation radius of a typical Be disc (20 i?*). 



locations of the electrons, seed photons and ob server can have an important 



effect, just as when computing 77 absorption ( Kirk et al. . 1999l h The high 



velocities of the shocked material (c/3 and higher in numerical simulations) 
can also cause Doppler (de)boosting effects. Finally, the pulsar wind itself 
is very likely anisotropic. 



3 Radio tails 

The radio emission in LS I +61°303 and LS 5039 has been resolved on 
angular scales of 10-100 mas and widely interpreted as du e to relativistic 



jet emission by analogy with compact jets in X-ray binaries (|Paredes et al 



2002; iMassi et all l2004h . The one-sidedness of the radio emission in LS 5039 



was interpreted as a Doppler effect, implying a jet speed of « c/3. Strong 
changes in the orientation of the (one-sided) radio emission in LS I +61°303 
on a timescale -C P or b were interpreted as jet precession. On larger scales, 
the radio emission appears to be more double-sided in LS 5039 and stable. 

Emission from cooled PWN particles can explain the resolved radio fluxes 
very well. Isolated ms pulsars interacting with the ISM are known to emit 
well-collimated cometary tails of cooling shocked material in the direction 
opposite to motion, extending to parsec scales. Furthermore, colliding winds 
in some Wolf-Rayet binaries produce non-thermal radio emissi on whose ap- 



pearance changes with orbital phase (D ougherty et a l.. 2003). In a 7-ray 
binary, the appearance of the PWN will combine the effects of cooling in a 
comet tail whose direction is changing with the pulsar's orbital motion. 

The aspect depends on orbital elements, shock geometry and cooling 
physics; but its basic features can be derived. For a face-on system (i=0°), 
the orbital motion creates a spiral of decreasing intensity with distance. The 
speed of the cometary material is typically crc so the step of the spiral is 
~ cc-Porb- When the spiral is seen edge-on (i=90°), the PWN has a double- 
sided, jet-like aspect with peak intensity alternating on both sides with P rb- 
Inclination effects are less pronounced in a highly eccentric system, which 
has a clear preferred direction for most of the orbit, with dramatic changes 
in the position angle of the tail when the pulsar moves around periastron 3 . 

The PWN model predicts that the radio morphology at high resolution 
(~ acP or b) should change on the orbital period. Maps at such resolutions 

3 This could explain the large P.A. changes seen bv lMassi et al.1 J2004) in LS I +61°303 
(e=0.72) compared to LS 5039 (e=0.35). Unfortunately, there are no high resolution maps 
of PSR B1259-63 (e=0.87) available for comparison. Both radio PWN should be easier to 
resolve than in LS 5039 as the orbital timescales are longer. 




Figure 3: Radio appearance of the PWN nebula of LS 5039 on milliarcsecond 
scales at </>=0. The basic shape is a spiral of decreasing intensity with a step 
crPorbC (~2 mas here with <r=0.01). The projection on the sky assumes i=60°. The 
curved jet-like emission reproduces very well the only available VLBI observation. 
The emission will move left and right with orbital phase as the pulsar follows its 
orbit. At resolutions ^> crP or bC, far-away particles from many orbits contribute, 
and the change in morphology is averaged out, leaving a steady radio source with 
an extension along the projected major axis. 



depend mostly on geometric orbital effects and less on details of cooling. The 
milliarcsecond LS 5039 map close to periastron is in excellent agreement with 
the published VLBI map (Fig. 3). At resolutions 3> acP ov \> the appearance 
is due to the summed emission of particles over many orbits (spiral steps), 
so the emission will be stable in morphology and intensity. Good knowledge 
of particle cooling and flow conditions are then required to make robust 
predictions at low resolutions. 



4 Conclusion 

Emission from a compact plerion provides a common, coherent framework to 
understand 7-ray binaries. These sources are the short-lived progenitors of 
accretion-powered HMXBs, which have slowly-rotating spundown pulsars. 
About 30 such sour ces are expected in our Galaxy given the present-day 



HMXB population (|Meurs & van den Heuvel ll989b 



High-energy emission occurs on small scales. Strong variations in the 
location of the pulsar wind termination shock (LS I +61°303, PSR B1259- 
63) and 77 absorption of VHE photons on starlight (LS 5039) can cause 
periodic changes at X-rays energies and above. Particles leave the vicinity 

Consistent with seeing 3 within 3 kpc 



of the pulsar as they cool, forming a cometary tail. Periodic changes in the 
aspect of this tail can be resolved in radio on scales acP ov ^, providing a test 
for the PWN scenario. The radio nebula is steady on larger scales. 

Compact plerions enable access to pulsar wind physics on hitherto in- 
accessible scales. Detailed understanding of the VHE emission, notably in 
LS 5039, can provide important clues as to the shock geometry, magnetic 
field intensity and acceleration processes around young rotation-powered 
pulsars. 
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